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Abstract 

Lie systems in Quantum Mechanics are studied from a geometric 
point of view. In particular, we develop methods to obtain time evo- 
lution operators of time-dependent Schrodinger equations of Lie type 
and we show how these methods explain certain ad hoc methods used 
in previous papers in order to obtain exact solutions. Finally, several 
instances of time-dependent quadratic Hamiltonian are solved. 

1 Introduction 

The use of tools of modern differential geometry has been shown to be very 
useful in many different problems in physics and in particular Lie groups and 
Lie algebras have played a prominent role in the development of Quantum 
Mechanics. The main concern of Lie was the integration of systems of dif- 
ferential equations admitting infinitesimal symmetries but as a byproduct of 
his work we have available a lot of relevant tools to deal with many different 
problems not only in differential geometry and classical mechanics but also 
in Quantum Mechanics. 
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Our aim here is to show the efficiency in solving quantum problems of 
the theory developed by Lie for dealing with systems of differential equa- 
tions admitting nonlinear superposition rules for solutions and that there- 
fore maybe considered as a generalization of non-autonomous linear systems. 
More specifically, we will be mainly interested in finding the time evolu- 
tion operator for quantum systems described by time-dependent quantum 
Hamiltonians which turn out to be the quantum counterpart of the above 
mentioned Lie systems of differential equations. 

Time-dependent Schrodinger equations of Lie type are Schrodinger equa- 
tions of motion for which the solutions of the equations can be obtained from 
the solutions of an equation in a Lie group G related to the Hamiltonian in a 
way to be explained in the paper. A particular example is the harmonic os- 
cillator, almost ubiquitous in physics, when the mass and angular frequency 
are not constant but become real time-dependent parameters p]-[5], may be 
as a consequence of interaction with the environment j6]. Methods for solv- 
ing such systems have been developed since long time ago, the generalized 
invariant method proposed by Lewis and Riesenfeld [7j being a typical and 
powerful method. Generalizations of such time-dependent harmonic oscilla- 
tors have also been considered and general quadratic systems have been very 
often studied from different perspectives [8] [T6j. 

This kind of equations appears very often in physics: quantum optics 
[T71 ITS] , quantum chemistry [T9], Paul trap [20]- [23], quantum dissipation 
[23], fluid dynamics [25] . etc. In particular, exact solutions of this type of 
equations appear many times in the literature [6]- [27]. The methods used to 
obtain the solutions of these problems are numerical or reduce the problem 
of finding out the solution of the time-dependent Schrodinger equation to 
solve certain differential equations using certain ad-hoc hypothesis. This is, 
for instance, the way used in the Lewis-Riesenfeld method [U [2H] or in the 
method of unitary transformations [Til 129] . 

Following the techniques of [30]- [35] we develop in this paper a geometric 
interpretation in which the solutions of this sort of Schrodinger equations 
are obtained in a natural way, by solving a system of differential equations 
like in other previous papers, but without any additional ad-hoc hypothesis. 
Our presentation is also an important improvement from the algorithmic 
point of view, it offers information about the difficulty of solving certain 
Hamiltonians, allows us to solve different physical examples at the same 
footing and provides us with a frame to explain different methods used in 
the literature. Finally, our method is an improvement because it shows that 
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Lie's ideas can be applied in the case of Schodinger equations. 

The paper is organized as follows: Section 2 gives a review of the mathe- 
matical theory of Lie systems and summarizes the result of the main theorem 
due to Lie [30] . Some simple examples of Lie systems are also given and spe- 
cial emphasis is made on the main property, the possibility of relating them 
with a particular type of equations on a group. Section 3 is devoted to 
present explicit formulas for a method of solving such equations which is a 
generalization of the method proposed by Wei and Norman for linear systems 
[36l [37] . The motion of a classical particle under the action of a linear poten- 
tial is analysed from this perspective. Lie systems in Quantum Mechanics are 
studied in Section 4 and the reduction method in Section 5 where the inter- 
action picture is revisited from a geometric point of view. Some applications 
of Lie systems in Quantum Mechanics are studied in Section 6, and more 
specifically, time-dependent quadratic Hamiltonians. Other examples of Lie 
systems appearing in the literature are pointed out in Section 7. Finally, the 
conclusions and outlook are presented in Section 8. 



2 Lie systems of differential equations 

In this section we will detail some known results about Lie systems that 
will be applied along this paper. First, we recall that time evolution of 
many physical systems is described by non-autonomous systems of differential 
equations 

— = X\x,t), i = l,...,n, (1) 

for instance, Hamilton equations, or Lagrange equations when transformed 
to the first order equations by considering momenta or velocities as new 
variables. In modern geometric terms, such a system is substituted by a 
t-dependent vector field 

n f> 

i=l 

whose integral curves satisfy ([I]). 

The theorem of existence and uniqueness of solution for such systems 
establishes that the initial conditions x(0) = (a; 1 (0), . . . , x n (0)) determine 
the future evolution. It is also well-known that for the simpler case of a 
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homogeneous linear system there is a (linear) function F : R n +n — > M n , 
given by 

X = F(x ( i),...,X( n ),/Ci,...,/c n ) = hx {1) H h k n X(n) , (2) 

in such a way that the general solution can be written as a linear combination 
of n independent particular solutions, X(i)(t), . . . , X( n )(t), 

X(t) = F(X(i)(t), . . . ,X( n )(t), h, ...,k n ) = kiX(i)(t) H h k n X(n){t) , (3) 

i.e. x(t) given by (j3J) is a solution for any choice of k = (kx, . . . , k n ) and 
for each set of initial conditions, the coefficients k = (ki,...,k n ) can be 
determined. In a similar way, for an inhomogeneous linear system, there is 
an affine superposition function F : R™(' 1+2 ) — > R" given by 

x = -FO(i), • • • ,X( n+ i), h, ...,k n ) 

= X(i) + h(x(2) - + . . . + k n (x [n+ i) - X(i)) , (4) 

and the general solution can be written as the corresponding affine function 
of (n + 1) independent particular solutions 

x(t) = F(x(i)(t),...,X( n+1 )(t),k 1 ,...,k n ) 

= X(i)(t) + h(x(2)(t) - X(i)(t)) + . . . + k n (x(n +1 )(t) - X(i)(t)) . (5) 

Under a non-linear change of coordinates both systems become non-linear 
ones. However, the fact that the general solution is expressible in terms of 
a set of particular solutions is maintained, now the superposition function 
being no longer linear or affine, respectively. 

The very existence of such examples of systems of differential equations 
admitting a non-linear superposition function suggested to Lie the analysis 
and characterization of such systems. He arrived in this way to the problem of 
characterizing the systems of differential equations for which a superposition 
function, allowing to express the general solution in terms of m particular 
solutions, does exist. The solution of this problem due to Lie [Stl] asserts 
that, under very general conditions, such systems are those which can be 
written as 

— = h(t)e i (x) + ■■■ + br{t)C{x) , i = 1, • • • , n , (6) 
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where b r functions depending only on t and £ m , a = 1, . . . , r, 

are functions of x = (x 1 , . . . , x n ), such that the r vector fields in M. n given by 

X a = £ ^(x 1 , .■-,*") ^ , a = 1, . . . , r, (7) 

i=i 

close on a real finite-dimensional Lie algebra, i.e. the vector fields X a are 
linearly independent and there exist r 3 real numbers, fa/3 1 , such that 

r 

[X a ,Xp] = ^/ a/3 7 X 7 . (8) 

7=1 

Even if the theorem was not stated with today level of rigour, the result 
is essentially true. For an intuitive geometric proof see [31], and for a more 
geometric approach and the uniqueness of such superposition rule see [3U [35] . 

From the geometric viewpoint, Lie systems are those corresponding to a 
t-dependent vector field which is a t-dependent linear combination of vector 
fields 



X(x,t)=Y l ba(t)X a ( 



a=l 

with vector fields X a closing on a finite-dimensional real Lie algebra. Many 
of the applications in physics and mathematics of such Lie systems have been 
developed by Winternitz and coworkers [38]-[46j. 

We have mentioned in this section that homogeneous linear systems are 
Lie systems. Let us consider such a system like 

^ = Y i A* j (t)xf t i = l,...,n, (9) 

i=i 

for which 

n f) 

x=Y, Ai >® x 'te> (10) 

which is a linear combination with time-dependent coefficients, 

n 

X=Y,A i j (t)X ij , (11) 
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of the n 2 vector fields 



• d 

X v =xJ g-i> *>J = 1> ■■■>"> ( 12 ) 



for which 

[Xij,Xki] 



3 9 I 9 



dx l ' dx k 



S a x i—-S kj x l — 
dx k dx l 



i.e. 



[X ij ,X kl \ = 8 u X kj -5^X il , (13) 

which means that the vector fields {X a = Xij,a = (i — l)n + j}, with 
i,j = 1, . . . ,n, appearing in the case of a homogeneous system, close on a 
n 2 -dimensional real Lie algebra isomorphic to the Ql(n, R) algebra. Actually 
they are the fundamental vector fields corresponding to the natural action of 
GL(n, M) on R n . 

Another remarkable example, with many applications in physics, is that 
of Riccati equation, which corresponds to n = 1 and m = 3 (see e.g. [HI H7J 

SB]). 

From the practical viewpoint, the most important property of Lie systems 
is that as the vector fields X a appearing in the t-dependent linear combina- 
tion defining the system 

r 

X = J2b a (t)X a , 

a=l 

are assumed to close on a finite-dimensional real Lie algebra g, if they are 
complete vector fields, they generate an effective action $ : G x R n — > R n 
of a connected Lie group G with Lie algebra jj, see jl9], on R n and if we 
determine a curve g(t) in G starting from the neutral element, g(0) = e, and 
such that 

r 

R 9-Ht)*g(t)9{t) = - ^2b a (t)a a = a(t) , (14) 

a=l 

where {ai, . . . , a r } is a basis of the Lie algebra g closing the same commuta- 
tion relations as the X a , then, the solution of ([6]) with initial condition x(0) 
is given by 

x(t) = ${g{t),x(0)) . 
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The equation flHj) is sometimes written with an abuse of notation as 



In this way, the problem of finding the general solution of ([6]) is reduced 
to that of determining the above mentioned curve in G solution of (jl~4"j) and 
starting from the neutral element [3TJ E2J ED]- In the particular case of 
linear systems the group is (a subgroup of) the general linear group GL(n, M.) 
and the action is linear, and therefore, <3>(g(t), ■) provides the time evolution 
operator. 

All the preceding process can be straightforwardly generalized to deal 
with Lie systems in a general manifold M, the obtained superposition rule 
being generally only local. 

The remarkable point is that once a Lie system with Lie group G in a 
manifold M which is a homogeneous space for G has been solved on its group 
G for a curve a(t) in g it is possible to obtain the solutions of all other Lie 
systems with the same Lie algebra g and curve a(t) in any other homogeneous 
space N for G. Moreover, if the Lie group involved in the problem is solvable 
the problem can be solved by quadratures independently of the curve given 
in g. Otherwise, the problem is to be solved, when possible, for each curve 
in q separately, and for some specific examples the solution may be explicitly 
shown. Of course, if a problem is solved for a certain Lie group G for all curves 
in g the problem is also solved for any problem given by a Lie subgroup H 



3 The Wei— Norman method 

In this section we will describe a method to solve directly the equation (|14[) 
which is a generalization of the one proposed by Wei and Norman [36j [37] for 
finding the time evolution operator for a linear systems of type dll(t)/dt = 
H(t)U(t), with £7(0) = I (see also |47j). We will only give here the recipe of 
how to proceed, the proof can be found for instance in [311 [32J, [50] • 

The generalization of Wei-Norman method consists on writing the solu- 
tion g(t) of (1141) in terms of its second kind canonical coordinates w.r.t. the 



V 




a=l 



of G. 



7 



basis {ai 



. , a r } of the Lie algebra q, for each value of t, i.e. 



r 



g{t) = Y[exp(-v a (t)a a ) = exp(-«i(t)o 1 ) • • ■ exp(-v r (t)a r ) , (15) 



a=l 



and transforming the equation (fl4l) into a system of differential equations for 
the unknown functions v a (t). The curve g(t) we are looking for is the one 
given by the solution of this last system determined by the initial conditions 
^«(0) = for all a = 1, . . . , r. The minus signs in the exponentials have been 
introduced for computational convenience. Now, it can be shown that using 
the expression fTTS]) and after some mathematical manipulations, equation 
( I14p becomes the fundamental expression of the Wei-Norman method [50] 



with v a (0) = 0, a = 1, . . . , r. The resulting system of differential equations 
for the functions v a (t) is integrable by quadratures if the Lie algebra is solv- 
able [36, 37], for instance, for nilpotent Lie algebras. Finally, this system of 
equations depends only on the structure constants of the Lie algebra. 

As an interesting example, from the physical point of view, illustrating the 
possible applications of the theory we can consider the motion of a classical 
particle under the action of a linear potential. Such example has been studied 
in [31] and is reproduced here for the sake of completeness. This model, with 
many applications in physics, has been often considered both in classical and 
quantum approaches (see, e.g. [511 [52]) and it has recently been studied by 
Guedes [28J; its solution using the theory of Lie systems was given in [50J. 
The classical Hamiltonian is 



For instance, when f(t) = qEo + qE cos cut, it describes the motion of a 
particle of electric charge q and mass m driven by a monochromatic electric 
field. 

The classical Hamilton equations of motion are 




(16) 



(17) 



V 

x = — , 
m 

V = -/( 



(18) 
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and therefore, the motion is obtained by two quadratures 

Pot 1 



x (t) = x + ^- - — / dt' / fit") dt" , (19) 
m m Jo Jo 

p(t) = po- f f(t')dt'. (20) 
Jo 

In the geometric approach, the t-dependent vector field describing the 
time evolution is 

m ox op 

which can be written as a linear combination X = — X± — f(t) X 2 , with 

V 9 Y 9 

Xl= P7T ' X2 = 7T ' 
ox op 

being two vector fields closing with X 3 = d/dx a three-dimensional real Lie 
algebra isomorphic to the Heisenberg algebra, namely, 

[X U X 2 ] = -X 3 , [X 1} X 3 ] = , [X 2 , X 3 ] = . (21) 

The flow of these vector fields is given, respectively, by 

0i (t, (x ,Po)) = (x +Pot,p ) , 
(p 2 (t, {x ,p )) = (x ,p + t) , 
<f) 3 (t, {x ,po)) = (x + t,p ) . 

In other words, this corresponds to the action of the Lie group of upper 
triangular 3x3 matrices, the Heisenberg group, on M 2 , 







1 1 




a 3 


(h 


\-\ 





1 








V o 





1 









u J 







Let {ai, ^3} be a basis of the Lie algebra with non-vanishing defining 
relations [ai,a 2 ] = —a 3 . Then, the corresponding equation in the group (fill) 
becomes in this case 

gg~ l = a x + f(t) a 2 . 



ni 
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Now, choosing the factorization g = exp(— u 3 a 3 ) exp(— w 2 a 2 ) exp(— U\ ai) 
and using the Wei-Norman formula fll6p we will arrive to the system of 
differential equations 

ui = — , u 2 = -f(t) , u 3 - Ui u 2 = , 
m 

together with the initial conditions 



«i(0) =w 2 (0)=« 3 (0) = , 

with solution 

Ul = - , U 2 = - f f(t') dt' , « 3 = -- rft' /' tZt* 

^ y m Jo Jo 

Therefore the motion will be given by 









H 









( xo \ 


)i 









1 i -iSl<XSU(t")dt" 

o i -Itf{t')dt> 

1 

which reproduces (|20I) . Inverting the matrix we will obtain two constants of 
the motion corresponding to the values of the initial conditions 



h = P(t)+ f f(t')dt' 
Jo 

h = x(t)--(p(t) 
m V 



f fit') dA t + - f df f fit") dt" , 

Jo J m Jo Jo 



the first one being the one given in [28] 



4 Schrodinger Lie systems in Quantum Me- 
chanics 

In this section we will review a way to generalise Lie's ideas in order to use 
them in the case of Schodinger equations [341 56J. As a new result we show 
how to apply our method to obtain solutions for quadratic Hamiltonians. 
Many particular cases of this Hamiltonian can be found in the literature, but 



10 



obtained by ad hoc or approximate methods as it will be explained in detail 
in a next section. In our method, nevertheless, we found an algorithmic way 
to find them all. 

As far as Quantum Mechanics is concerned, let us first remark that 
the separable complex Hilbert space of states TL can be seen as a (infinite- 
dimensional) real manifold admitting a global chart [22] ■ The Abelian trans- 
lation group provides us with an identification of the tangent space T^TL at 
any point G TC with TC itself, the isomorphism being obtained by associating 
with ip ETC the vector tp G T^TC given by 

i>f (<!>)•= (jtf^ + t^) > 

for any / G C°°{H). 

Through the identification of TL with T^TC at any <fi G TL a continuous 
vector field is just a continuous map A: TC —>■ TC. In particular, a linear 
operator A on K is a special kind of vector field. Usually, operators in 
Quantum mechanics are neither continuous nor defined on the whole space 

n. 

The most relevant case is when A is a skew-self-adjoint operator, A = 
—i H. The reason is that TC can be endowed with a natural (strongly) sym- 
plectic structure, and then such skew-self-adjoint operators are singularized 
as the linear vector fields that are Hamiltonian. The integral curves of such 
a Hamiltonian vector field A = —i H are the solutions of the corresponding 
Schrodinger equation [53]. Even when A is not bounded, it can be shown 
that if A is skew-self-adjoint it must be densely defined and its integral curves 
are strongly continuous and defined in all TC. 

On the other side, the skew-self-adjoint operators considered as vector 
fields are fundamental vector fields relative to the usual action of the unitary 
group U {TC) on the Hilbert space TC. 

It is important to remark that given a r-dimensional real Lie algebra of 
skew-self-adjoint operators —iH a 

[iH a , iHp] = c aP 7 iH 7 , c a p 7 G M (22) 

we can choose a basis {a a \ a — 1, . . . , r} of an abstract Lie algebra $j isomor- 
phic to that of the X a such that the Lie brackets of the elements a a of this 
Lie algebra, denoted by [•, •] satisfy 

[a a , a p ] = c aP 7 a 7 , , c a/3 7 G JR. (23) 
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Now, Lie system theory applies to the case in which a t-dependent Hamil- 
tonian can be written as a linear combination with t-dependent real coeffi- 
cients of some self-adjoint operators, 



H(t) = J2 b a(t)H a , (24) 

cc=l 

where the Hamiltonians H a are such that the skew-self-adjoint operators 
—iH a close a real finite dimensional Lie algebra under the commutator bracket 
as indicated in (1221) . 

When this happens, we can associate a Lie algebra q acting on 7i by 
skew-self-adjoint operators in such a way that there exists a basis {a a } of q 
with —iH a being the fundamental vector field associated with a a . 

The linear map X : g — > X(7i), X : a a i— > —iH a is a Lie algebra isomor- 
phism. Then, when is given by (1241) . the Schrodinger equation 

^ = -iif = - h *{t)iH a il> (25) 



can be reduced to an equation on G like (T141) determined by the curve in 

r 

given by a(i) = — ~S^b a (t)a a . Once equation (TT4T) with initial condition 



a=l 



p(0) = e has been solved we can obtain the general solution of the Schrodinger 
equation. 

As a first instance, we can see how the method works in the very simple 
case of a quantum time- dependent linear potential, which has recently been 
studied by Guedes |28j, instead of using the Lewis and Riesenfeld invariant 
method [7]. Our method is an improvement with respect to previous ones, 
because it allows us to obtain the solution in an algorithmic way. 

In this case the quantum Hamiltonian is 



H 9 = ^ + f(t)X. (26) 

In this quantum problem, as pointed out in [MlIM], the quantum Hamil- 
tonian H q may be written as a sum 

H q =-H x -f(t)H 2 , 
m 
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with 

H\ = — , Ho = — x . 

2 

But —i H\ and —i H2 are skew-self-adjoint and close on a four- dimensional 
Lie algebra with — i H 3 = —i P, and — iH^ — i I, which is an extension of the 
opposite of the Heisenberg Lie algebra (121]) . i.e. the defining commutation 
relations are 

[-i Hi, "« H 2 ] = -iH 3 , [-1 H u -1 H 3 ]=0, [-1 H 2 , -1 H 3 ] =-iH±. 

Therefore, the Lie algebra of Schrodinger equation given by the Hamil- 
tonian H q is like in the analogous classical Lie system. Despite of the fact 
that H q given by (j26|) is time-dependent, it is a Lie system and thus we 
can find the time-evolution operator by solving a related equation on the 
corresponding Lie group by the Wei-Norman method. 

More explicitly, let {ai, a 2 , a 3 , a^} be a basis of the Lie algebra with non- 
vanishing defining relations [01,02] = —03 and [02,03] = —04. The equation 
f fT4|) in the group to be considered is now 

Rg-u g 9 = -— 01 + fit) o 2 . 
m 

In order to find the expression of the wave-function in a simpler way, it 
is advantageous to use the factorization 

g = exp(— 1> 4 a4)exp(— 1> 2 02) exp(— 1> 3 a 3 ) exp(— v\ a\) . 

In such the Wei-Norman method provides the system 

vi = — , V2 = ~f{t) , 
m 

1 1 2 

v 3 = — v 2 , v 4 = -—v 2 , 
m 2m 

that, jointly with the initial conditions fi(0) = f2(0) = ^3(0) = f4(0) = 
determines the solution 

vi(t) = ^, v 2 (t) = - J\t' f(t') , (27) 
v 3 (t) = -- f dt' f dff(f), (28) 



^0 



2m 



-t' - 2 



v 4 (t) = - J- / dt'l I dt"f(t") ) . (29) 



\Jo 
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Then, applying the evolution operator onto the initial wave-function ijj(p, 0), 
which is assumed to be written in momentum representation, we have 

Mp) = U(t,0W (p) 

= exp(it> 4 (t))exp(iv 2 (t)X)exp(— ivs(t)P)exp(— ivi(t) P 2 / 2)%jj (p) 

= exp(^ 4 (t))exp(^ 2 (t)X)e~ i(,,sWp+,;iWp2/2 Vo(p) 

= e -i(-V4(t)+v 3 (t)(p-v 2 (t))+v 1 (t)(p-v 2 (t)) 2 /2)^ p _ V2 (+X\ ) 

where the functions Vi(t) are given by (127|) . fl28l) and (1291) . respectively. 

We can proceed in a similar way with the quadratic Hamiltonian in the 
quantum case, given by [55] (see [34"] ) 

H q = a(t) ^ + (3{t) XP + PX + 7 (t) + 5 (t) P + c(f ) X + (30) 

where X and P are the position and momentum operators satisfying the 
commutation relation 

[X,P] = iI. 

It is important to solve this quantum quadratic Hamiltonian because of its 
appearance in many branches of physics. 

The Hamiltonian can be written as a sum with t-dependent coefficients 

H = a(t) Pi + (3(t) P 2 + 7 (t) P 3 - 6{t) P 4 + e(t) P 5 + <f>(t) P 6 

of the Hamiltonians 

P 1 (x,p) = ^, H 2 (x,p) = ^(XP + PX), H 3 (x,p) = ^, 

P 4 (x,p) = -P, H 5 (x,p)=X, H 6 (x,p)=I, 
which satisfy the commutation relations 





iH 2 \ 


= %H X , 


[iH 2 ,iH 3 ] 


= iH 3: 


[iH 3 , 


iH 4 ] 


= iH b , 


'iH x , 


iH 3 ] 


= 2iH 2 , 


[iH 2 ,iH A ] 






iH 5 ] 


= 0, 


iH 1: 


iH 4 ] 


= 0, 


[iH 2 ,iH 5 ] 


= ~H 5 , 








iH x , 


iH 5 ) 


= -iH A , 













and [iH a , iH 6 ] = for a = 1, . . . , 5. 
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This means that the skew-self-adjoint operators —i H a generate a six- 
dimensional real Lie algebra. Thus, we can define a basis {ax, . . . , a^} with 
the same structure constants as the iH a with Lie product the commutator 
of operators. 

This six-dimensional real Lie algebra is a central extension of the Lie 
algebra arising in the classical case by the one-dimensional Lie subalgebra 
generated by a 6 . It is a semidirect sum of the Lie subalgebra sl(2, R) spanned 
by {ai, a 2 , 0,3} and the Heisenberg-Weyl Lie algebra generated by {04, 05, a^}, 
which is an ideal. 

In full similarity with the classical case, in order to find the time-evolution 
provided by the Hamiltonian (I3"U|) we should find the curve g(t) in G such 
that 

6 

R g -i* g g = - ^ b a (t) a a , flr(O) = e, 
«=i 

with 

61 (t) = a(t) , 62 (t) = Pit) , 6 3 (t) = 7(*) , h{t) = Sit) , b^t) = e(t) , 6 6 (t) = <f>(t) . 

This can be carried out by using the generalized Wei-Norman method, 
i.e. by writing g(t) in terms of a set of second class canonical coordinates. 
For instance, 

g(t) = exp(-f 4 (t)a 4 )exp(-t;5(t)a5)exp(-t; 6 (t)a 6 )x 

xexp(— vi(t)a 1 )exp(—v 2 (t)a2)exp(— v 3 (t)a3) (31) 

and a straightforward application of the above mentioned Wei-Norman method 
technique leads to the system 

v 1 = b 1 + b 2 vi + b 3 v\ , v 4 = 64 + - b 2 v A + b x v 5 , 
< i> 2 = b 2 + 2 63 vi , i> 5 = 65 - 63 t> 4 - \ b 2 v 5 , (32) 



1 

2 



with ^(0) = U2(0) = v 3 (0) = u 4 (0) = u s (0) = v 6 (0) = 0. 
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If we consider the following vector fields 



9 9 1 2 9 

X 2 = Vi— + — + -W 4 - o^5^— , 

afi 9t> 2 2 9t> 4 2 9t> 5 

2 9 9 d 9 1 n 9 

X 3 = ^7— + 2 Wl — + e^ 2 — - y 4 — + ^l—, 
av\ 0V2 9f 3 9f 5 2 9t> 6 



i 



9t> 4 ' 

9 9 



(33) 



X 5 = a ^4 

9^5 OVq 

X - 9 
9f fi 



it is a straightforward computation to see that these vector fields close on 
the same commutation relations as the corresponding {a a } and thus (1321) is 
a Lie system related with the same Lie group as the Hamiltonian (1301) or its 
corresponding equation in its Lie group. 

Now, once the v a (t) have been determined, the time-evolution of any state 
will be given by 

\ip{t)) = exp(-v A (t)iH 4 )exp(-v 5 (t)iH 5 )ex.-p(-v 6 (t)iH 6 )x 

x exp(-%(t)iifi)exp(-f; 2 (t)i/J 2 )exp(-'y 3 (t)ii?3)|V'(0)) . 



and thus 

\if)(t)) = exp(v 4 (t)iP)exp(—v 5 (t)iX)exp(—v e (t)iI)x 

( P 2 \ ( PX + XP\ ( X 2 \ 
x exp I -v x {t)i— J exp ( -v 2 (t)i J exp I -v 3 (t)i— J \ijj(Q)) , 

5 The reduction method in Quantum Mechan- 
ics 

We start this section with a quick review of the reduction technique explained 
for example in [321 SB] and then we obtain some new results. First, as a new 
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improvement, while in some before works like [3"2l US] some sufficient con- 
ditions where explained to perform a reduction process, here we show that 
these conditions can be considered as necessary. Also, we will use the reduc- 
tion technique to explain the interaction picture used in Quantum Mechanics 
and we will review from the point of view of our theory the method of unitary 
transformations. 

It has been proved in Section 2 that the study of Lie systems in homoge- 
neous spaces can be reduced to that of the solution of equations 

r 

Rg-U g 9 = b a (t)a a = a(t) E T e G (34) 

a=l 

with g(0) = e. 

The reduction method developed in [32J has also shown that given a 
solution x(t) of a Lie system in a homogeneous space G/H, the solution of 
the Lie system in the group G, and therefore the general solution in the given 
homogeneous space, can be reduced to that of a Lie system in the subgroup 
H. More specifically, if the curve g(t) in G is such that x(t) = <&(g(t),x(0)) 
with $ being the given action of G in the homogeneous space, then git) = 
g(t)g'(t), where g'(t) turns out to be a curve in H which is a solution of a Lie 
equation in the Lie subgroup H of G. Actually, once the curve g(t) in G has 
been fixed, the curve g'(t), which takes values in H, satisfies the equation 

This transformation law can be understood in the language of connec- 
tions. It has been shown in [321 [56] that Lie systems can be related with 
connections in a bundle and that the group of curves in G, which is the 
group of automorphisms of the principal bundle Gxl [5B] , acts on the left 
on the set of Lie systems in G, and defines an induced action on the set of 
Lie systems in each homogeneous space for G. More specifically, if x(t) is 
a solution of a Lie system in a homogeneous space M defined by the curve 
a(t) in q, then for each curve g(t) in G such that g(0) = e we see that 
x'(t) = $(g(t),x(t)) is a solution of the Lie system defined by the curve 

a'(t) = Rg-i, s g + Ad(g)a(t) . (36) 
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which is the transformation law for a connection. 

In summary, the aim of the reduction method is to find an automorphism 
g(t) such that the right-hand side in fl36l) belongs to T e H = f) for a certain Lie 
subgroup H of G. In this way, the papers [32j [56] gave a sufficient condition 
for obtaining this result. In this section we shall study the above geometrical 
development in Quantum Mechanics and we find out a necessary condition 
for the right-hand side in (13"6"|) to belong to f). 

Lie systems in Quantum Mechanics are those such that 

r 

H(t) = ^2b a (t)H a . (37) 

o=l 

with —iH a closing under the commutator on a finite dimensional real Lie 
algebra D. Therefore, by regarding these operators as fundamental vector 
fields of an action of a connected Lie group G with Lie algebra g isomorphic 
to o, we can relate the Schrodinger equation with a differential equation in 

r 

G determined by curves in T e G given by a(t) = —'^^b a (t)a a by considering 

a=l 

—iH a as fundamental vector fields of the basis of q given by {a a }. 

Now, the preceding methods allow us to transform the problem into a 
new one in the same group G, for each choice of the curve g(t) but with a 
new curve a'(t). The action of G on 7i is given by an unitary representation 
U, and therefore the time-dependent vector field determined by the original 
Hamiltonian H will become a new one with Hamiltonian H' . Its integral 
curves are the solutions of the equation 

^ = -iH'W (38) 

where 

- iH'{t) = -iU{g{t))H{t)U\g{t)) + U (g(t))U\g(t)) (39) 

That is, from the geometric point of view we have related a Lie system 
on the Lie group G with certain curve a(t) in T e G and the corresponding 
system in 7i determined by a unitary representation of G with another one 
with different curve a'(t) in T e G and its associated one in 7i. 

Let us choose a basis of T e G given by {c a \ a = 1, . . . , r} with r = dimjj, 
such that {c a \ a = 1, . . . , s} be a basis of T e H, where s = dim f), and denote 
{c a | a = 1, . . . , r} the dual basis of {c a \ a = 1, . . . ,r}. In order to find g 
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such that the right hand term of (1361) belongs to T e H for all t, the condition 
for g is 

c a (Ad(g)a(t) + Rg-i^g) = , a — s + 1, . . . , r . 

Now, if 8 a is the left invariant 1-form on G induced from c a the previous 
equation implies 



djr 1 

dt 



( R- g -i* e a(t) - ) = , a = s + l,...,r. (40) 



Let g = g" 1 , the last expression implies that Rg* e a{t) — g is generated 
by left invariant vector fields on G from the elements of F). Then, given 
7r L : G — ► G/H, the kernel of 71"^ is spanned by the left invariant vector fields 
on G generated by the elements of f). Then it follows 

4(i^ e a(t)-|) = 0. (41) 

Therefore, if we use that 71^ o = —X^ o 7r L , where X^ denotes the 
fundamental vector field of the action of G in G/H and X^ denotes the 
right-invariant vector field in G whose value in e is a a , we can prove that 
7r L {g) is a solution on G/H of the equation 

^=twy(j))- (42) 

Cf=l 

Thus we obtain that given a certain solution in f) related to the 
initial g(t) by means of according to g(t) = g(t)g'(t), then the projection 
to G/H of g{t) is a solution of (|42|). 

The result so obtained shows that whenever ^'(t) is a curve in if then 
g{t) verifies (14"2|) . Moreover, as it has been shown in [32], if ^(t) satisfies (T4"2|) . 
then (yf'(t) is a curve in H satisfying (135^) . The previous result shows that 
such a condition for obtaining fl35l) is not only sufficient but necessary too. 
Thus, we obtain a new result which completes the one found in [32J. 

Finally, it is worthy of remark that even when this last proof has been 
developed for Quantum Mechanics, it can be applied also in ordinary differ- 
ential equations because it appears as a consequence of the group structure 
of Lie systems which is the same for both the Schrodinger and general dif- 
ferential equations. 
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5.1 Interaction picture and Lie systems 

As a first new application of the reduction method for Lie systems, we will 
analyse in this section how this theory can be applied to explain the inter- 
action picture used in Quantum Mechanics. This picture has been shown to 
be very effective in the developments of perturbation methods. It plays a 
role when the Hamiltonian can be written as a sum of a simpler Hamiltonian 
Hi and another perturbation Hi. In the framework of Lie systems we can 
analyze what happens when the Hamiltonian is 

r r 

H = Hx + V(t) = H x + J2 b *(t)H a = h(t) = 1, (43) 

a=2 a=l 

where the set of skew-self-adjoint operators {—iH a | a — 1, . . . , r} is closed 
under commutation and generates a finite dimensional real Lie algebra. The 
situation is very much similar to the case of control systems with a drift term 
(here Hi) which are linear in the controls. The functions b a (t) correspond to 
the control functions. 

According to the theory of Lie systems, let us consider a basis {a a | a = 
1, . . . , r} of the Lie algebra with corresponding associated fundamental vector 
fields —iH a . The equation to be studied in T e G is ( 1341) and if we define 
g'(t) = g(t)g(t), where g{t) is a previously chosen curve, it holds a similar 
equation for g'(t) given by (13H . 

If in particular we choose g(t) = exp(ait) we find the new equation in 
T e G 

Rg>-i* g <g' = -Ad(exp(ait)) ^2b a (t)a^j = -exp(ad(ai)t) ^^& Q (t)a Q j . 

(44) 

Correspondingly, the action of G on 7i by a unitary representation de- 
fines a transformation of Ti, in which the state ip(t) transforms into ip'(t) = 
exp(iHit)il)(t) and its dynamical evolution is given by the vector field corre- 
sponding to the right hand side of (jHJ). In particular, if {02, . . . ,a r } span an 
ideal of the Lie algebra q the problem reduces to the corresponding normal 
subgroup in G. 
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5.2 The method of unitary transformations 

A second application of the theory of Lie systems in Quantum Mechanics and 
in particular of the reduction method is to obtain information about how to 
proceed to solve a Lie Hamiltonian in Quantum mechanics. Even when the 
theoretical results used here can already be found in the literature, the way 
to proceed can be considered a new improvement. 

Every Schrodinger equation of Lie type is determined by a Lie algebra 
g, a unitary representation of its connected Lie group G in 7i and a curve 
a(t) in T e G. Depending on g there are two cases. If g is solvable we can use 
the reduction method in Quantum Mechanics to obtain the general solution. 
If g is not solvable, it is not possible to integrate the problem in terms of 
quadratures in the most general case. But it may be possible to solve a 
problem completely for some specific curves as for instance it happens for 
the Caldirola-Kanai Hamiltonian [24J. A way of dealing with such systems 
is to try to change the curve a(t) into another one a'(t), easier to handle, as it 
has been done in the previous subsection of the interaction picture. In a more 
general case than the interaction picture, although any two curves a(t) and 
a'(t) are always connected by an automorphism, nevertheless, the equation 
which determines the transformation can be as difficult to be solved as the 
initial problem. Because of this, it is interesting to look for a curve easy 
enough to be solved but that we can connect to the initial problem. In any 
case, we can always express the solution of the initial problem in terms of a 
solution of the equation determining the transformation. In certain cases, for 
an appropriate choice of the curve g(t) the new curve a'(t) belongs to T e H for 
all t, where if is a solvable Lie subgroup of G. In this case we can reduce the 
problem from g to a certain solvable Lie subalgebra f) of g. Of course, in order 
to do this, a solution of the equation of reduction is needed, but once this is 
known we can solve the problem completely in terms of it. Other methods 
have alternatively been used in the literature, like the Lewis- Riesenfeld (LR) 
method. However, this method seems to offer a complete solution only if g is 
solvable. If g is not solvable, the LR method offers a solution which depends 
on a solution of differential equations, like in the method of reduction. 

To sum up, given a Lie system with Lie algebra g, with G acting unitarily 
on TC and determined by a curve a(t) in T e G, the systematic procedure to be 
used is the following: 

• If g is solvable we can solve the problem easily by quadratures as it 
appears in E2] 
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• if g is not solvable, we can try to solve the problem for a given curve 
like in the Caldirola-Kanai Hamiltonian in [21], by choosing a curve 
g(t) transforming the curve a(t) into another easier to solve, like in 
the interaction picture. If this does not work we can try to reduce 
the problem, like in the time-dependent mass and frequency harmonic 
oscillator or quadratic one-dimensional Hamiltonian in [U I2S1 ESI EZ] 
to an integrable case. 



6 Applications in Quantum Mechanics 

In this section we shall apply our methods to obtain time-dependent evo- 
lution operators of several problems found in the physics literature in an 
algorithmic way. In particular, we will analyze several examples of quadratic 
Hamiltonians. These examples are studied in the literature under different 
approaches but here we will study them using the same viewpoint. In our 
description we will classify the developed examples in terms of whether the 
related Lie group is either solvable or not. Also, in the non-solvable cases, 
we will describe some approaches to the study of these Lie systems. 



6.1 Solvable Hamiltonians 

Quadratic Hamiltonians describe a very large class of physical models. Some- 
times, one of these physical models is described by a certain family of quadratic 
Hamiltonians that can be considered as a Lie system related with a Lie sub- 
group of the one given for general quadratic Hamiltonians. When this Lie 
subgroup is solvable the differential equations related with it through the 
Wei-Norman methods are solvable too and the time-evolution operator can 
be explicitly obtained. In this subsection we will deal with some instances of 
this case. In this cases, we can find the explicit solution of these problems 
already the literature by using for each case a different method. 

First, we will fix our attention at the motion of a particle with a time- 
dependent mass under the action of a time-dependent linear potential term. 
The Hamiltonian that describe this physical case is 

H -^T) +S(t)x - (45) 

The Lie algebra associated with this example is a central extension of the 
Heisenberg Lie algebra and at the same time, a Lie subalgebra of the one 
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obtained for quadratic Hamiltonians. A basis for the Lie algebra of vector 
fields related with this physical model is 



p2 

Z x = i—, Z 2 = iP, Z 3 = iX, ,Zi = iI (46) 

which closes on a Lie algebra under operator commutation, 

[Z x , Z 2 ] = 0, [Z u Z 3 ] = 2Z 2 , [Z u Z A ] = 0, 

[Z 2 ,Z 3 ] = Z 4 , [Z 2 ,Z 4 ]=0, (47) 

[Z 3 , Z A ] = 0. 

This Lie algebra is solvable, and then, the related equations obtained through 
the Wei-Norman method, can be solved by quadratures for any pair of time- 
dependent coefficients m{t) and S(t). The solution of the associated Wei- 
Norman system allows us to obtain the time-evolution operator and the wave 
function solution of the time-dependent Schrodinger equation. 

This Hamiltonian has been studied in [58J for some particular cases using 
ad-hoc methods and in general in [32]. Here, we will use the Wei-Norman 
method. As a quadratic Hamiltonian its equation in the group G is a par- 
ticular case of the one related with (I30I) 



1 

m(t) 



R g -u g g = — 177IT 0,1 ~ = a>Ms(t) ■ (48) 



If we use the factorization given in fl3TT) 

g(t) = exp(-v A (t)a4)exp(-v 5 (t)a 5 )exp(~VQ(t)a 6 )x 

x exp(— vi(t)a 1 )exp(—V2(t)a 2 )exp(—v s (t)a s ) (49) 

we can solve the Schrodinger equation by the Wei-Norman method through 
the set of differential equations (I3"2|) for this particular case 

1 

vi = —77^ , v 2 = 0, v 3 = , 
m(t) 



2 



v 4 = ——, v 5 = S(t), v e = -S(t)vi 



m(t) ' w ' ' w 2m(t) ' 

with initial condition t>i(0) = v 2 (0) = v 3 (0) = v^O) = v 5 (0) = v 6 (0) = 0. 
The solution of this system can be expressed using quadratures because the 
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related group is solvable: 



v 2 (t) = 
v 3 (t) = ^ 

^ = f^tu) ( f S{v)dV 
v 5 (t) = J S{u)du 

vM - -i' sM (Sw) U? (a)dw )) du - /'^ ( rs(v,di - 2 

(50) 

and the time-evolution operator is obtained by using the last expressions in 



U (g(t)) = exp(— t> 4 (t)2if 4 )exp(— t> 5 (t)iiJ 5 )exp(— t> 6 (t)ii/ 6 )exp(— Vi{t)iHi) 
= exp(t> 4 (t)iP)exp(— t> 5 (t)zX)exp(— f 6 (t)i/)exp(— iv\(t)P 2 /2) 

(51) 

Now, in the case of constant mass, (150{) reads as 
t 

Vl (t = - 
m 

v 2 (t) = 
v 3 (t) = 

f 4 (t) = — / I / S(v)dv j du 



v§(t) = I S{u)du 



«*(*) = (^)/ {l S ^ dw ) dv ) du -^l {I S{v)dv ) du - 

which gives the time-evolution operator if we use them in fl5T|) 

Now, as we have obtained the time-evolution operator for a Hamiltonian 
related with any curve ciMs{t) we can consider particular instances of it. For 
example, for the curves with constant mass m and S(t) = qe + qe cos(cot) 
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studied in BE! we obtain 



v 1 (t) = —, v 2 (t) = Q, v 3 (t) = 0, 
m 

vAt) = -(2e + e uj 2 t 2 — 2ecos(u;t)) , vAt) = —(e ujt + esm(u>t)) , 

2muj z uj 

and 

-a 2 

v 6 {t) = ~7T)~ 3 (4e^ 3 t 3 - 3e(e - Ae )ut + 

3e(4e + 2e (aA 2 - 2) - 3e cos(cut)) sin(cut)) . (52) 

The way to obtain a solution with arbitrary non-constant mass and S(t) = 
qe + qecos(ut) was pointed out in [23] and solved in [52] • From our point 
of view, the most general solution is straightforward from (15 Up because all 
cases in the literature are particular instances of our approach with general 
functions m(t) and S(t). 

Now, we can obtain the wave function solution of this system. We know 
that the wave function solution with initial condition ip is 

1h(x) = U(g(t))4> {x) 

= exp(iv 6 (t))e^(-vAt)iP)exp(-v 5 (t)iX)exp ^--v^i^j ip (x) 

(53) 

However, expressing the initial wave function ip in the momentum space as 
4>o(p) the solution is given in the same way as before but with U(g(t)) in the 
momentum representation. In this case the solution with initial condition 

Mp) is 
Up) = u( g (t))Mp) 

= exp(-w 6 (t))exp(w 4 (t)iP)exp(-?;5(t)iX)exp (-^(t)^) Mp) 

( P 2 \ 

= exp(— zf 6 (t))exp(f 4 (t)iP)exp(— f 5 (t)zX)exp I —ivi(t)— (fio(p) 



/ (p 4- vAt))" 2 
exp(-ive(t))exp(v 4 (t)iP)exp I -ivAt) - ) <fio(p + v 5 (t)) 

exp -iv 6 (t) + ivAt)p - ivt(t) ) Mp + M*)) 
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6.2 Non-solvable Hamiltonians and particular instances 

In subsection 16 . 1 1 the differential equations associated to the Quantum Hamil- 
tonians treated there were Lie systems related with a solvable Lie algebra. 
Thus, by the Wei-Norman method, it has been shown that they were inte- 
grable by quadratures. When this does not happen it is not easy to obtain a 
general solution for the differential equations. Now, we will describe some ex- 
amples of this type of quadratic Hamiltonians. In general we will not obtain 
a general solution in terms of the time- dependent functions of the quadratic 
Hamiltonians. Nevertheless, we will show that for some instances of them 
the differential equations can be integrated. Note that explicit solutions of 
these Hamiltonians cannot generally be obtained as our theory explains, but 
under some integrability conditions on the coefficients the solution 

can be worked out. 

As a first case consider the Hamiltonian of a forced harmonic oscillator 
with time-dependent mass and frequency given by 

This case, either with or without time-dependent frequency, has been 
studied in [U [281 EH] ■ The equations which describe the solutions of this Lie 
system by the method of Wei-Norman are 

«! = —— + m (t)u 2 (t)vl 
m(t) 

i)2 = 2m(t)tu 2 (t)v i 
v 3 = e V2 m{t)uj 2 (t) 
1 

m[t) 

va = fit) - m(t)uj 2 (t)v4 

v 6 = Kn{t)u 2 {t)vl - f(t)v 4 - 7^^ v l 

with initial conditions fi(0) = f2(0) = ^3(0) = ^4(0) = fs(0) = vq(0) = 0, 
where we have used the factorization fl3Tl) . The solution of this system cannot 
be obtained by quadratures in the general case because the associated Lie 
group is not solvable. Nevertheless, we can consider a particular instance of 
this kind of Hamiltonian, the so called Caldirola-Kanai Hamiltonian [241. In 
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this case, for the particular time-dependence m(t) = e rt m , u>(t) = u and 
f(t) = 0, the Hamiltonian is given by 

H = ^e rt + lm e~ r W X 2 . (54) 

In this case the solution is completely known and is given by 

2e rt 



m (r + lu coth (^o)) 
t> 2 (t) = rt + 2 logu; — 2 log sinh ^7^0^ + w cosh 



2 



r + u coth(^tuo) 



^2' 

Vi{t) = 0, v s (t) = 0, v 6 (t) = 0, 

where a>o = \J r 2 — 4c<Jq. This example shows that the problem may also be 
solved exactly for particular instances of curves in q of Lie systems with non 
solvable Lie algebras. Another example is the following one 

P 2 muj\ 
2m + 2{t + k) 2 ' 
for which the solution of the Wei-Norman system reads 

2(k + t)((k + t)*° - fc°°) 



H=— + nu ■ " „ X 2 , (55) 



Vl< ^ ~ m{k^{u Q -l) + {k + t)^{u)Q + l))' 
v 2 (t) = (l+uj )log(k + t) - (l + LJo)\ogk + 2\og(2k 0o u ) 

-2 log(r°(o) - 1) + (k + t)°°{u + 1)) , 

2mu 2 (k + t)" - k^ 

V *"' ~ k k°o(u Q -l) + (k + t)^(u + l) 
v 4 (t) = 0, v 5 (t) = 0, v 6 (t) = 0, 



where now cu = a/1 — Aoj 2 . 

Other examples of Hamiltonians, which can be studied by our method, 
can be found in [ 24] . We just mention two examples which can be completely 
solved 

P 2 1 

Hi = — + -m {U + Vcos{u t))X 2 , 
2mo 2 

H 2 = ^e rt + ±m e- r W X 2 + f(t)X. 
27 



The first one corresponds to a Paul trap as has been studied in J5U], and 
admits a solution in terms of Mathieu's functions. The second one is a 
damped Caldirola-Kanai Hamiltonian analysed in [58] . 



6.3 Reduction in Quantum Mechanics 

Quite often, when a non-solvable Lie algebra is involved in a quantum prob- 
lem, it is interesting to solve it in terms of (unknown) solutions of differential 
equations. Next, we study some examples of how to proceed with the method 
of reduction in order to deal with problems in this way. So, we will obtain 
that the reduction method not only can be applied in Quantum Mechanics 
but also allows us to solve certain problems in an algorithmic way. As far as 
we know, this kind of application of the theory of reduction in Lie systems 
to a Schodinger equation is new. 

Consider an harmonic oscillator with time-dependent frequency whose 
Hamiltonian is given by 

H=^ + ^(t)X 2 . (56) 

As a particular case of the Hamiltonian described in Section H] this example 
is related with an equation in the connected Lie group associated to the 
semidirect sum of st(2,K) with the Heisenberg Lie algebra generated by the 
ideal {a^, 05, as} 

R g -u g g = -01 - Q 2 (t)a 3 , g(0) = e. (57) 

Since the solution of this equation starts from the identity and {ai,a 2 ,a 3 } 
close on a s[(2,R) Lie algebra, then the Hamiltonian H of (156]) is related 
with the group SL(2,R). Actually this is due to the isomorphism of such 
a group with the symplectic one, and in higher dimension the group will be 
Sp{2n, R) instead of SL(2n, R). 

As a particular application of the reduction technique we will perform the 
reduction from G = SL(2, R) to the Lie group related with the Lie subalgebra 
(3 = (ai). To obtain such a reduction, we have shown in Section [5] that we 
have to solve an equation in G/H, namely 

^ = EMW(J)) (58) 
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where are the fundamental vector fields of the action A of G on G/H. 
Now, we are going to describe this equation in a set of local coordinates. 
First, in an open neighborhood U of e £ G we can write in a unique way any 
element of SL(2, M) as 



g = exp(a 3 a3)exp(a 2 a 2 )exp(aiai). (59) 



where we choose 



ai= (o J)' a2 = Ko -i)' a3= (-°i o)' (60) 

This decomposition allows us to establish a local diffeomorphism be- 
tween an open neighborhood V C Gj H and the set of matrices given by 
exp(a3a3)exp(a 2 a2). Now, the decomposition (1591) reads in matrix terms as 

a(3\_fio\fe o \ / i i) \ _ f e o \ / 1 v 
7 5 ) ~ \(f) i;v° ^ / V 1 / V w ^ / v 1 

If we express <f), 8, ip in terms of a, (3, 7 and 5 we obtain 



(61) 



a P\_( 1 \ f a \ f 1 (3/a \ _ f a \ ( 1 (3/a 
7 5 J ~ \ -f/a 1 J \ a' 1 J \ 1 J ~ \ 7 a' 1 J \ 1 

(62) 

Thus, we can consider the projection tt l : U C G — > G/if given by 



_L 



7T 



a /3 \ [a 



7 5 / 17a 1 



F (63) 



which allows to represent locally the elements of G/H as the 2x2 lower 
triangular matrices with determinant equal to one . Now, given \ g : g'H £ 
G/H — > gg'H £ G/H as \ g o^ L = tc l oL g then the fundamental vector fields 
defined in G/H by a\ and a 3 through the action of G on G/H are given by 



t=0 



,Mexp(- tal )(^))^(7 7 ; a2 



7r L I exp(— ta 3 ) 



a /3 \\ _ / 
7 5 J J — 1 ct 



t=o 

and the equation in V G G/H is described by 

a \ _ / -7 
7 — da -2 y I f2 2 (t)a 7a" 



(64) 



(65) 
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and thus we want to obtain a solution of the system 

y 1 (66) 

7 = —a . 

Then, if ot\ is a solution of the system (166]) the curve git) that verifies g{t) = 
g(t)h(t), where hit) is a solution of an equation defined in the group related 
with f) = (ai), reads 

m = ( a i _-. ) (e^) 



and the curve which acts on the initial equation in SL(2,M.) to transform it 
into one in the mentioned Lie subalgebra is given by g(t) = 

g{t) = exp(— 2 log «i a 2 )exp ( — — a 3 ) = exp(— a 1 d 1 a 3 )exp(— 2 logaii a 2 ) 

V «i / 

(68) 

and this curve transforms the initial equation in the group given by (1571) into 
the new one given by (cf. fl36l) ) 

a 'W = 277\ a i ' ( 69 ) 

which corresponds to the Hamiltonian 

H' = - P 2 (70) 

and the induced transformation in the Hilbert space 7i that transforms the 
initial H into if' is 

exp (i^PX + IP)) exp • (71) 



Both results can be found in [57] . 

There are other possibilities of choosing different Lie subalgebras of g in 
order to perform the reduction, however the results are always given in terms 
of a solution of a differential equation. 
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7 Conclusions and outlook 



It has been shown that the geometric approach to Lie systems of differ- 
ential equations can be extended to the framework of quantum mechanics 
and several examples of application have been developed. In particular we 
have studied time-dependent quadratic Hamiltonians for a single particle and 
other related models. Some of them are involved in models for dissipative 
Quantum Mechanics or harmonic oscillators in external fields. In all of these 
applications we have developed methods to obtain exact solutions that sum- 
marize many of the different techniques used in the literature in a simple and 
unifying framework. The method of reduction has also been revisited and 
a new theoretical result and some new applications in Quantum Mechanics 
have been obtained. 

But these are not the only applications of the methods developed here, 
they can be used in many other applications. For example, Hamiltonians of 
the following type are studied in Quantum Optics [HI] 

H = uoia\ai + u] 2 a\a,2 + flo(t)(e lujt a\a 2 + e~ lujt a\ai) 

+ n L (t)a\ + fi£(t)a a + Q R (t)at + n* R (t)a 2 . 

In this case, it can be shown that this Hamiltonian is a Lie system and thus 
we may study systems of interacting harmonic oscillators in external fields in 
our formalism. Many particular cases of this Hamiltonian have been studied 
in several papers and in different contexts. 

Another different example is a system of coupled spins in magnetic fields 
like the following one found in [S3] 

H = -2A Si • Sj + B(t) • S i 

i,j i 

In this example our method allows us to obtain exact solutions when B(i) 
is a three dimensional vector field in the form 

B(t) = (Bsin6cos(ut),Bsm9sm(ut),Bcos{9)). 

On the other hand, our viewpoint allows us to study important properties 
of these and other physical systems too. For example, the last example 
is related with Berry's phases. Our method allows us to obtain the time- 
evolution operator. Then this information can be used to determine Berry's 
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phases. More specifically, the integrability conditions obtained in [25] can be 
generalised to solve exactly some Hamiltonians for certain time-dependent 
coefficients. We can also try to apply control theory to determine the relation 
of Berry's phases with the field B(t), study the adiabatic approximation, etc. 
All these topics are important and will be analyzed in forthcoming papers. 
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